


‘ 


nd LIPA 






Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1971 


Identification for linear electrical power 
system models 


DeVille, Thomas G 


: Massachusetts Institute of Technology 


http://ndl.handle.net/10945/15746 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
(8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist | | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


1 i KNOX appointed — and published -- scholarly author. 

ih : Dudley Knox Libra Naval Postgraduate Schoo 
LIBRARY dley b ry i | g d hool 

http://www.nps.edu/library 






411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


IDENTIFICATION FOR LINEAR ELECTRICAL POWER 
SYSTEM MODELS 


Thomas G DeVille 








(sei CeT ION FOR LINSAR 
EIECTRICAL POWER SYSTEM MOREIS 
by 
THOMAS GEORGE DEV TLLS 
B.S., United States Coast Guard Academy 


(1966) 


SUeieesD IN PARTIAL BUMPILMEIT 
OF THE REQUIREMeW?S FOR THE DEGREFS OF 
MECHANICAL ENGITEER AND 
MESTER OF SCIENE 
at the 
PESShCE Cael IS INSTITUTE OF TECHNOLOGY 


May, 1971 





LIBRARY 
NAVAL POSTGRADUATE SCHOOL 


MONTEREY, CALIF. 93940 
=Oe 


IDEN IFICATION FOR LIMEAR 


RIFCTRICAL POWER SYSTEM MODELS 


by 
THOMAS GEORGE DEVILIE 


Submitted to the Department of Mechanical Engineering and the 
Department of Electrical Engineering on May 14, 1971, in partial 
fulfillment of the requirements for the cegrees of Mechanical Engineer 
anc Mester of Science in Electrical Engineering. 


ABSTRACT 


Analyses of electric power systems ordinarily include only a 
small portion of the entire network. Effects external to the network 
of interest are modeled as loads or generators. A method is proposed 
for identifying an equivalent model of the external network without 
taking any measurements in the external network itself. The method is 
simple and could easily te implemented to supplement state estimation 
techniques. 


A linear load flow model which relates real power to voltege phase 
angle is Cceveloped for enalyzing the entire network and from this two 
models are Ceveloped for relating measurements within the system of 
interest to the parameters of the equivalent model. It is found that 
there are two components to the equivalent model. One is a fixed struc- 
ture and the other is an equivalent power. Since both models for 
identifying the equivalent network are in the form of an input, output 
relation with an unknown additive disturbance, a general solution for 
identifying linear static systems is found subject to certain conditions. 
Applying the solution to the two models for identifying the equivalent 
network, it is found that one of the mocels hes an inherent tendency to 
produce consistently biased estimates of the equivalent model perameters. 
However, that seme model hes the acvantege that it uses data for its 
input, output that is readily available arc relatively accurate. 


A simulation made to test the proposed method is discussed for a 
Sample system and the results agree well in both the actual paremeters 
identified and predicted error. Several important points are discussed 
with a view toward implementation. 
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1. INRCODUCTION 


Analyses of electrical power transmission networks usually carry 
the implication that the network under consideration can be isolated 
from all other networks. In actuality this can seldom be realized 
since the network is interconnected with one or more outside networks, 
and whatever happens to the network being analyzed is determined at 
least to some extent on how it interacts with the others. Therefore, 
it would be convenient if the effects of all the outside networks on 
the network of interest could be summarized by one equivalent network. 
In most cases the equivalent network would be considerably smaller than ~ 
all the actual outside networks that it would represent. For instance, 
if the New Englard network were interconnected with the rest of the 
United States at a half dozen points, es far as the New England network 
would be concerned, the entire rest of the United States could be 
shrunk to an equivalent network having a half dozen buses. A method is 
proposed for finding the equivalent network which could be used for 


contingency planning and state estimation techniques [3 Su, [s Jj, 


fsl,ft9}, and fre]. 
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1.1 General Problem: Identificetion of Linear, Static Svstems 

When analyzing a system, the usual approach is to surround the 
system of interest with a conceptual boundary. That within the boundary 
is idealized in some respect to emphasize certain features. Further, 
the idealized system is allowed to communicate with the outside only 
through certain conceptual ports. This idealized representation of 
actuality is called the system model and by analyzing the model conclu- 
sions can be drawn about the actual system. If the contents of the 
idealized system and what enters through the ports are completely known, 
then what happens to the idealized system can be found at least in 
principle. If the model is static, then only what currently enters is 
necessary to know the current state of the svstem. In such a case, if 
only knowledge of the current state is required, then the effect of the 
system outside the bourdary is summarized by what enters through the 
ports. However, to know what will happen in the future, ae necessary 
to know what will enter in the future, and it may be that what enters 
through the ports is at least partly determined by how the system inter- 
acts with the outside. If everything outside the system is also ideal- 
ized, then everything of interest can be modeled by the two idealized sys- 
tems. The first, the original system model, represents what is within 
the boundary. The second, the external system mcdel, represents all 
that is not represented by the first. Further, the two models can only 
interact through the originally specified ports. The external system 
model may be rather extensive and it might be desirable to find an 
equivalent external model which would interact with the original system 


model in the same way. 
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The process of finding the equivalent system model, called identi- 
fication here, is imvestigated for a class of linear, static system 
models. More specifically, it is assumed that the model consists of a 
structure which remains constant and two types of variables, inputs (or 
causes), and outputs (or effects). The mathematical notation using 


matrices is simply y = C 


u where u represents an input vector, y repre- 
sents an output vector, and © represents the structure. 

If the original system model with its conceptual ports to the 
outside is placed in this framework, then it might be possible to 
Separate what exists at or passes through the ports into the same two 
types of variables. For the model of a physical system, the product of 
the two variables frequently represents a generalized power passing 
through the ports. An exemple would be the voltage level existing at 
an electrical terminal and the current pessing through the terminal. 

By combining inputs from the ports with inputs internal to ohe system, 
the outputs can be evaluated. In the same way knowledge of inputs at 

the ports, inputs within the external system model, and the structure of 
the external model allow outputs of the external model to be fourd. That 
is, each model may be analyzed iniepemently of the other by knowirg 

what passes through the ports. However, in order to predict future 
outputs of the original system, knowledge of how the two models interact 
is necessary. That is the purpose of identifying the equivalent external 
system model. 

The inputs and outputs can be divided into those belonging to the 
original system, those of the external system, ard those common to both. 
Two general aporoaches to identifying the equivalent model are investi- 


gated. One involves combining both the original and external models 
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and then identifying the new model using all the inputs not in the 


external model and only the outputs common to both. Another approach 
is to use the inputs and outputs at the ports or terminals to identify 
the equivalent external system separately. Both approaches are analo- 
gous to a black box which contains both the external model structure 
and inputs belonging to the external mcdel. Inputs from the original 
system enter the black soem outputs are those which are common to 


beth 
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1.2 Specific Problem: Electrical Power Transmission Networks 


Ordinarily when analyzing electrical power transmission systems, an 
arbitrary boundary is drawn about some part of the entire network and 
that part of the system is studied assuming power flowing in through 
the trarsmission lines on the boundary is known. The specific problem 
investigated is that of determining how an interconnected electrical 
power system affects a mer eine subsystem. The subsystem will be 
nemed "Own System" and referred to as (© for convenience. The rest will 
be named "External System" and referred to as XS. Together (CS and XS 
form the entire system named "Whole System” or WS. The External System 
may be the entire outside world as viewed from one region. Or CG may | 
be a higher voltage transmission network looking down into a more com- 
plex lower voltage distribution network. XS may also be a relatively 
extersive region where power and voltage measurements are not made which 
is surrourded by 0S where measurements are made. To aid in analysis 
notation CS is further subdivided into "Internal System” or BD which 
has no immediate connections to 4S and "Boundary System" which has 
immediate connections to XS. What is to be determined is an equivalent 
model for XS named "Equivalent System" or ES which will affect CS the 
same as XS, 

It will be assumed that measurements are available in © but not from 
XS. There are two aspects which set this particular identification prob- 
lem apart from many others. One is that the identification must be 
primarily passive. It will not be possible to make any great manipula- 
tion of input signals (power flows) that unduly upset the system or do 


not conform to power consumption requirements. Also it can be expected 
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eee 
that power requirements in XS will be about the same as in (%S, and 
since the unknown inputs in XS will cause unknown changes in the 


outputs of OG, "signal to noise ratios" will be about one to one. 
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1.3 Contents of Thesis 

The basic conventional model which is used for analyzing networks 
will be developed in Section 2. Variables of interest for this model 
are real power, reactive power, voltage magnitude, and the relative 
voltage phase angle. The model is nonlinear and involves cross coupling 
of all variables. However it will be shown that under certain conditions 
the effect of voltage magnitude changes on real power is much less than 
the effect due to phase angle changes, so that an approximate linear 
model can be developed which relates real power to voltege phase angle. 
The linear model is not conventional but it is not new either. One of 
its chief advantages is a very significant reduction of computation 
requirements. But even though a linear model will ce used for identi- 
fying the model for Equivalent System, it is still possible to incor- 
porate the identified equivalent linear model in the nonlinear model for 
Cwn System. In Section 3 the linear model for External eyenen is 
reduced to an equivalent model. Then two models for use in identifying 
Equivalent System are ceveloped using the equivalent model. One of the 
models is attractive from the point of view that it only requires obser- 
MametOors aw the tie lines joining G to AS, but it will turn out that 
there are inherent problems with this model due to the correlation 
among variebles at the tie lines. Both models will be placed in a 
general form and in Section 4 the solution to the general problem will 
be found sudject to certain conditions. The solution, which has a very 
Simple form, is derived from both an assumed mathematical probability 
model for the disturbances and the method of weighted least squares. 


Error analysis equations are also derived to study the source of errors 





ale a 
and find a measure of confidence in the identified parameters. 
Applying the solution of the general problem to the two specific 
identification models in Section 5, the effects of the conditions 
placed on the solution are analyzed. One shortcoming of the more 
easily implemented model is discussed in terms of independent inputs. 
£ simulation made to test the proposed methods of identification is 
discussed in Section 6, and the results show excellent correspondence 
with the theory. Errors in the identified parameters are within pre- 
dicted accuracy. Finally, in Section 7 several aspects of the specific 
problem are discussed which involve increasing identification accuracy 
and verification of the identified model. t is also pointed out how 
models for identifying the equivalent system can be valuable by them- 


selves for the purpose of network model reduction. 
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2. BiECTRICAL POWER TRAISMISS ION MOEIS 


2.1 Nonlinear load Flow Model 

In reality an electrical power system consists of a large variety 
of devices, mamy capable of storing energy with time constants and 
natural frequencies spanning a very wide range. However for the pur- 
pose of computing power flows in sinusoidal steady state only a few 
types of components need be considered. The most important of these are 
the trarsmission lines which transmit power frem points of generation 
to loads, the transformers, and the buses which form connecting points 
or nodes for transmission lines, transformers, generating stations, am 
loads. The values of interest are the power flows and bus volteges. 
Although the transformers, transmission lines, loads, and sources are 
capable of storing energy, when the system is in sinusoidal steady state, 
the usual approach is to treat energy storege devices as complex | 
impecances or edmittances and the sources am loads as imjectiors of 
complex power. Hence although the system is oscillating, it can be 
analyzed as a static network. 

A transmission line is of course a cistributed parameter device. 
However es with most such devices, it is possible to model it as a set 
of lumped elements provided the physical size is small compered to the 
wavelength of the power transmitted. Since the wavelength of 60 hertz 
alternating current is roughly the distance between New York and San 
Francisco, while most transmission lines are far less, lumped element 
models are quite valid for this purpose. The pi model, common for this 


type of cistriduted parameter device, is used here arm! consists of a 


resistance and inductance in series with shunt capecitance at each em. 





a7 = 





Figure 2.1 Pi Transmission Line Model 


eT YY Ts 
+O | = 


Figure 2.2 Electrical Power Network Lumped Model 
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Since the system studied is assumed to be in sinusoidal steady 
state, the capecitences and inductance can be represented as admit- 


tances. A bar above a symbol will denote a complex quantity (e.g. 
i jD 
Ey =a e =E, cos D, + JE, sinD,). An underlined symbol will 
denote a matrix or vector. Iet 
et Reo aeics _ 1 
ik ik Rey + Job) 


be the complex admittance of the series resistance and inductance, and 


let 
a, eee 
YS = IPS 


be the admittance of the shunt capacitance. The current entering the 


c ee: : th eee a ea = = ee 
transmission line at the i--° node is IL, = E, ys,, + (E, - Ey) Vans 

a ee = 
Complex power is P+ jQ =E I (I is the complex conjugate of I) so 
the power entering the transmission line from the ple node is 


ae SS Siav= =e 
PLiy, ~ SQL, = E, Ty, = 2; By ys; + Ey (EF; - By) yy, 


(2.1) 
The network model consists of admittances which represent the 

transmission lines and transformers connected to buses or nodes. Also 

connected to each node is a current source which represents the load 

or generator at the bus. Power injections into the buses are positive 

for generators and negative for loads. &s shown in Appendix A, a bus 

admittance matrix can be formed which represents the relation between 


the current sources at each bus to the bus volteges. Iet E 


ws pone 


vector whose elements are the complex bus voltages, se: a vector of 


the corresporming values of the complex bus current sources, and Be 





= 
the complex bus admittance matrix. Then the relation is 


N 
— i 
Lee ec Sus or” for’ tne “1 bus , I, e Yi; where there are 


N buses. Corresporndingly, the complex power injected into the i 


bus 
by the current source at the bus is 
N -jD. N jD - je 
: — Te i k ik 

Pemion ~E. > Y.,E, = 5,e eres ac 
az i iS ikk i Kel k ik 

N -j(6,,4+D.-D,) 

ae a Rit (2.2) 
wa) 2 k ik 


The load flow problem is that of given bus power injections, firm 
the bus voltages am the transmission line power flows. It should be 
noted that even though the network is modeled with linear elements, the 


power flow equations are nonlinear. 
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2.2 Linear Real Power - Voltage Angle Model 
Under certain assumptions it is possible to reduce the complexity 
of the model. From Equation (2.1) line real power flow is 


A lee -J APP) t 


ie 
PLey Real 1 YSsy + Ey Y se - E Ese 


_y 2 @ 
= Es Yay COS By ~ BiB, C08 (85, + Dy - DY) 


PL,, = E.y,,cos Bsr \Es - EL cos (D, - D,)|+ EEVy,,Sin 85, sin(D, - Dy) 
(2.3) 
If it can be assumed that voltage megnitudes do not vary much fron 
nominal operating points, that the transmission line reactance is much 
greater than the resistance so that Oy ~ 90°, and that voltage angle 


differences (D, = DL) are small, then when E, = (E) wominal? and 


BE, = (25) ominal? PL. a ES oY 3,5in, - Dy). For ee ante 
will be assumed that Ey = E. = 1.0. Using the approximation for the 
sine, sin(D, = D,) =D, -D, for D, - D, small 

PL, * vay, =o) (2c) 


As shown in Apperdix B, for such a case a linear load flow equa- 


bb 
tion, a = Loses? relates real power injections to voltage angles. 
For a system with (K+1) buses, Pi, is a Kxl vector of bus real power 


nie ; ; y : 
injections, D. . is a KxlI vector of bus voltage angles, and tous iS 2 


KxK matrix whose elements are nearly the same as the magnitude of 


elements of the complex matrix Y, The Ga? bus is the reference 


wns 
bus at which the voltage angle is specified. Power injection at the 


reference bus is such that the algebraic sum of real power injections 


of all buses is zero. Rom [6 3 and Rom and Schveppe {8} proposed the 
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use of this linear load flow model for on-line estimation of voltage 
angles and transmission line power flows while Baughman and Schweppe [2 J 


showed how its results compare with the nonlineer model. 


2.3 Discussion of the Tinear Load Flow Model 

Although the linear load flow model does not give the exact values 
of voltage angles obtainable by the nonlineer model nor dees it account 
for voltege magnitude variations, there ere some great advantages to 
using it for real power flow analyses even when the conditions on low 
line resistance and small voltage magnitude spread are only mildly met. 
Foremost is the simplicity of calculation. For (K+1l) buses the non- 
lineer model requires an iterative solution of a set of 2K simultaneous 
ronliresr equations. The linear model involves K simultaneous linear 
equations. If the Sewton-Rephson method is used to solve the first, 
some comparison between the two can be made. The Newton-Raphson 
method for one iteration is of the form Aa” = B” ac”. Since calcula- 
tions for solutions of simultaneous linear equatiors are about propor- 
tional to the square of the number of equations, one iteration for the 
nonlinear model requires about four times the calculations for the 
linear model. Also the Be matrix must be calculated at each iteration 
whereas daa in the linear model remairs fixed for a given system. 
Finally, cepending on the accuracy required, the Newton-Raphson method 
requires about four or more iterations. Another acvantage is the ease 
with which the linesr equations can be manipulated. This will become 
very important in transforming the network for finding an equivalent 


system as it results in much insight. Historically, before digital 
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computers were in wide use, load flow analyses were made with calculating 
boards which used physical components to model the network. The more 
versatile A-C calculating boards could account for voltage megnitudes 

and line power losses but were expensive anc there were only about 50 

in the United States O11). The linear load flow model is somewhat 
analogous to the D-C calculating boards which were more numerous amd 

less expensive. The nulieataas load flow model eccounts for individual bus 
real power injections while some D-C calculating boards used only one 
source. 

The error involved in using the linear load flow model of course 
depends on how well the three conditions are met. The assumption of 
small voltage angle differences between the ends of a transmission line 
usually results in the least error. Except in rather rare circumstances, 
angle differences are within 30° and usually they ere within 10°. Com- 


pering the values of the sine of an angle and the angle: 


Argle (degrees) Angle (radians) Sine ~ 
T 
10° 1327 Orn 745 Of 1736 
T 
20° 9. = 0.349 O32 
T 
30e pa Cee 0. 5CO 


Even up to 30° the approximation is within 5%, and under 10° it is 
Wwaroan 0.5%. 
The error involved in the other two assumptions can be analyzed 


through the nonlinear equation for line power flow. If the line 


conductance, Bs, = Ys, C08 Boys and the line suseptamrce, “bay =e sin Bers 


are used then equation (2.3) is 





a. 
PL, = 2,85,(F, - §,, cos, - D,)| + E,E,b,, sin(D, - D,) (2.5) 


In orcer to make a rather rough estimate of errors let "a" represent 


either Bey Or Ysy° Then when PL.» is approximated by b,,(0; - Dy) or 
¥ 5,0, - Dy) for small angle differences 

Pp ~ - se 

PL. B85, (5, E 1.) + E gEyb ics DL) 


bed the relative errer fom Pl... # Dae mis 
ik a K 


Ree 
relative error = PL ee 


2s ~ Py) - Byeay (By - By) - EVEyDy, 2, - Dy) 


ieee Oe (D, =o) 


ik ik 
— a BLED, 1 gy, By - Ey 
. : Dp. =p 
eb eee © ak 
Hor E. a 1.0 ard Ey ~ 1.0, € is approximately 
; Saeed 1 te 
; ae i Be, - 
ik 7 ieee) eck 


If the voltage drop is expressed as a percentage, AE = 1CO(E . = Ey )s 
and if angle cifference is expressed in degrees, AD = 57.3(D, - Dy) 


~ 50(D, - Dy), the relative error is approximately 


a of eammee 310 
= eee 
= ( £5, 2 »., AD aa 
ik sek 
The purpose of equation (2.6) is to show how the greatest portion of the 


error enters the linear trarsmission line power flow equation (2.4). 


The term a - BEY shows how voltege megnitude changes from nominal 
ik 





eA = 
1 3x AE 
values affect the error, while = =" = shows the effect of resis- 
2 Day AD 
tance em voltage drop. Table 2.1 compares line flows as computed 
both by the nonlinear and linear equations along with actual errors 
and approximate errors predicted by equation (2.6). It is not immedi- 
ately obvious whether Bey Or Ysy (actually BED Or EE ay where 
‘ z : = 2 
E, and E) are nominal values) should be used for a in PL. a(D, D,)- 
In view of the relatively small error and much less involved 


calculations required, and to preserve insight, the linear load flow 


model is used for the purpose of identifying the equivalent system model. 
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Se eee TON TOE I ICATION OF THE EQUIVAIE NT SYSTEM 


3.1 Tie Line Power Flow Model 

fs previously mentioned, two approsches to finding the equivalent 
system model will be used. The tie line power flow model views the 
external system as a separate system and attempts to find an equivalent 
model by using the tie line power flows and bus voltage angles where the 
boundary arourd (© crosses the tie lines connecting CS to XS. The bus 
injections and corresponding bus voltage angles can be separated into 
3 groups, those belonging to IS (Internal System), those belonging to 
BS (Boundary System), and those belonging to XS (External System). 
Collectively IS and BS form CS (wn System). PI and DI are the bus 
real power injection vector and the bus voltage angle vector respectively 
whose elements belong to IS. Similarly PB and DE belong to BS, and Px 
and DX belong to XS. The letters I, B, and X will usually refer 
values belonging to I, BS, and XS Peer elaa since by earlier 
definition I and XS are not immediately connected to each other by 


transmission lines, the linear load flow equations can be written: 





PI aig Ge 0 ||) Dr 
PB | = | XIB' YBB yYBx | | DB (3.1) 
P Deen | | bx 

(The transpose of a mstrix will be denoted by Y' = transpose of Y). 


The zero submatrices in the corners are cue to the lack of direct 
coupling between IS and KS. The only coupling between CS ard XS is 
reflected by the submatrix YBB. Since the elements of YEB are values 


of admittances of transmission lines which have at least one of their 








a 
ends connected to buses in BS, YBB can be separated into that which 
Betenes to Oogg2ei>, andethat which belongs to XS, YTL, so that 
YBB = YB@ + YIL. YL will be a diagonal matrix whose elemerts will be 


the admittances of the tie lines connecting G to XS. Equation (3.1) 


can be written 


PI i 1B 0) ep eB 9. oF | Bt 
Pay] = | ym" yecs oO} | DB] + {|O YrL yBx| | DB 
0 OQ oO} | Dx QO YBX' yxx| | Dx 
or 
PI 1D IB |e) DI 0 
= + (3.2) 
PB YIB' YBOS| | DB YTL DB + YBX DX 
PX = YBX'DB + YXX DX (3.3) 
let PPL = YTL DB + YBX DX (3.4) 


Tie line power flow, PIL, is a vector of length equal to P3 and repre- 
sents the power which flows out of the buses belonging to BS into the 
transmission lines which connect BS to XS. Solving for DX in equation 


(2.3) and substituting into equation (352) 


me = Soe “ema xx | 


7 
a 4 


PX 


*yBX'DB + YXX7" Px) 


i 


PEL = YEE DB + YSX(-YXx" 


(YL - YBX YXx"7YBX')D3 + YBX YXX""PX 


a 
2 —|— ——— 


- 


Define an equivalent external system bus admittance matrix, YEQX, and an 


equivalent external bus power injection vector, FEQX, such that 
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YEOK = -YBX YXXx77YBX! (3.5) 
FEQX = AQX Px (3.6) 
Agx = -YEx x”? (3.7) 


Then equations (2.2) and (3.4) can be rewritten 


PI YII yB] [DI 

= (3.8) 
PB - PIL YB‘ yBos! | DB 
PIL = (YTL + YEQK)DB - PEQX (3.9) 


Equation (3.8) is simply the linear load flow model for OS. By including 
PTL with FB the problem can be solved es usual. Equation (3.9) for PIL 
shows two contributions. The first, (YIL + YEQX)DB is the part of PIL 
Cue solely to the values of voltage angles existing at the buses of BS. 
The second, PEQX is the equivalent power injections in XS as seen by 


BS. It is the bus power injections from PX that are routed to BS by the 


Serucrure of %>. It would Ge desiraple to know both YEQX “ard PES but 





even knowledge of just YEQX would be valuable. For example, if the 
voltage angles and line power flows in CS were wanted after a change in 
bus power injections in CS, then if PX (and thus PEOX) did not change 


between times t, and t,, the change could be found by substituting 


iT oN 
equation (3.9) into (2.8) so that at any time t 


PI(t) YII YB DI(t) 


Pett) eee) | BY «(YBBR + YROX)| | DB(t) 


and for the change with PEQK (t..) = PROX (t, ) 





PI(t,)- PI(t,) Mit YI 


IS 


(t,) - DI(t,) 


)- BB(t,) YIB' (YBB + YEQX)| | DB(t,) - DB(t,) 


where DI(t,) and pB(t,) ere found by measuring FIL(t, ) and using 
equation (3.8). On the other hand, if the voltage angles and line 
power flows in © were wanted after a change in the status of one or 


more transmission lines in 6, then if YI, vB, and YBS are the admit- 


tance matrices after the change, for PI(t,) PI(t,), FB (t,) = Pa(t,); 





and FEQX (t) = PEQK (t, ) 


Yit YB DI(t,) YIT YB DI(t,) 
YIB' YBB+YEQX pB(t, ) YIB' YBB+YEQX DB(t,) 


Everything is known except BI(t,) and ba (t,) since egain DI(t, ) and 
DB (t,) can be found by measuring PIL(t,) and using equation (3.8). 


As formulated, in equation (2.9) DB is an input vector and PIL is 





an output vector while PEOX night be thought of as an unknown distur- 
barce. If a series of input, output measurements were taken, then it 
might be possible to find YEQX. Equation (2.9) will be referred to as 


the tie line power flow model. 
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3.2 F2oundsry Bus Impedance Model 
The second approech involves an impedance matrix instead of an 
admittance matrix. In order to obtain the required form, an equivalent 
internal system as seen by BS looking into IS will be fourd which is 
similar to the equivalent external system. The top row of equation (3.8) 
is 


Sie =vit Dr # vip DB 


Solving for DI 


= .V “1 


B+ yr 


ie 
4 
im 
3 
Ss 
) 
td 
KH 
a 


and substituting into the second row of equation (3.8) 


PB - PTL = YIB'(-YIT--YB DB + YIT PI) + YE 


i 
< 
le 


Using equation (3.9) for PIL and combining terrs 


+ PEOX - yinvir ie CS + YTL + YEQX - 117 yp )bB (3,10) 


Vv 


3 
Ld 


Define an equivalent internal acmittance matrix, YEQI, and an equi- 





valent internal bus power injection vector, PEQI, similer to YEQX am 


FRQX by 
YEQI = -YByir YB (3.11) 
PEQI = AQT PI (3 ale) 
HOI = -YIB¢vIT (G56) 


ome a ee 


Then equation (3.10) can be rewritten. 


PEQI + PB + PEQX = (YEQI + YBB + YEQX)DB (2205) 





= |= 


td 


QI is the power injection of PI which is routed to the buses of BS 


| 


by the strueture of I5. YEQI is the structure of IS as seen by BS. 

YBB was previously separated into that belonging to OS, YBOS and that 
belonging to XS, YTL. If YBCS is further separated into that belonging 
Womloyey ots ara thay belonging to BS, YBSS so that YBOS = YBIS + YBBS 


or YBB = YEIS + YBRS + YTL, then the quantity 
(YBIS + YEQI)D3 - PEQI 


is the power flowing out of the buses of BS into the transmission lines 
connecting IS and BS. If BS had no transmission lines unique to itself, 
then YBBS would be zero. Ey rewriting equation (3.14) with 


YBB = YPIS + YBBS + YTL 


PEQI + PB + PEQX = (YBIS + YEQI)DB + YBRS DB + (YTL + YEQX)DB (3. 14a) 


See ee 





Equation (3.14a) can be considered to be two matrix Thevenin or Norton 
Squivellemt circuits, one for 5 amd one for XS, coupled to the circuit 
for BS as illustrated in Figure 3.1. If equation (2.14) is solved for 


DB, then 


D3 = (YEQI + YBB + YEQX)~!(PEQI + PB + PEQX) 





Poeeeoe regedit 8, Galo + Hee + YEOOa de that part of the Whole 
System (WS) impedence matrix which relates PB to DB. Therefore denoting 


ZB8B for that part of the impedance matrix 
BB = ZB3(PEQI + PB) + ZBB PEQK (3.15) 


In this formuletion (PEQI + PB) is an input vector, DB is an output 


vector, and Z2B PEOX is an unknown disturbance vector. Again if a 








PEQI + PB + PEQX = (YEQI + YBIS)DB + YBES DB-+(YEQX + YTL)DB_ 





Mueure o.l Equivalemt Network Viewed From BS 


PT aa aaa 






ro 


TL = (YEQX + YTL)DB - PEOX 








Pigure 5.2 Tite“Line Power rlow Model 








—— = mem 


DB= ZEBB(FEQI + PB) + Z2BS FEQX 





Mweure 5.5 Beundary Bus Impedance Model 
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series of input, output measurements were teken, it might vce rossible 
moO findyeee. 

For comparison purposes and in order to be abdle to investigate 
sources of approximation error later, the analogous equations for the 


nonlinear model are Ceveloped in Appendix C. 


3.3 Separating Cwn System Estimation From Equivalent System Identification 


The tie line power flow model and the bourdary bus impedance model 
are iilustrated in Figures 3.2 and 3.3. As shown in Tadle 3.1, both 
approaches result in equations of the form z=Hu+v where u is the 
input vector, z the output vector, v a vector of unknown disturbarces, 
and H the structure matrix to be identified. Both z ard u were assumed 
to be known. However, since the state of CS (i.e. the bus voltages) can 
be estimated by knowing what power enters from XS through the tie lines, 
the problem of estimating the state of OS can be separeted from identi- 
fying the parameters of the Demet system. What happens to © 
depends on XS, but the effect of XS is summarized by the power flowing 
through the tie lines. Hence existing state estimation techniques (5 ], 
[8], and [9 Jean be applied to obtain z and u. In practice there will 
be some error essociated with metering power and estimating voltage 
angles. In such a case let We be the error of output estimate and Me 


be the error of input estimate. Then 


i) eae) 


- 


or Peewee (weaey + Uv) 
ee ee 2 te 





Out put 
Input 
Noise 


Structure 
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Table 3.1 


Identification Models 


General 
Model 


ic ic In 


i 


Pourdery Bus 
Impedance 
Model 


DB 


(P2+PEQT) 


ZBB PEQK 


Z3B 


Tie Line 
Power Flow 
Model 
PLL 
DB 
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By replacing v-v +8 


. rw 
we en with a new unknown disturbance vector v 
As 
Bs 


z=Hut 


the same approaches can be used. The only difference is that there is 


generally more disturbance associated with identifying H. 


oe Use of Changes in Variables 

Up to now z, u, and v have not been associated with time. From 
now on z(t), u(t), and v(t) will be used to designate a set of 2, u, 
and v for a particular time t. The disturbance vector v(t) will be 
treated as a noise term since u(t) is analcgous to an input Signal ard 
z(t) is analogous to an output signal corrupted by noise. Ordinarily 
when working with estimation or identification problems, it is much 
easier if noise terms have a zero mean, In the two approaches oe 
neither will in general have zero temporal mean noise rector vie). «Tt 
is expected that v(t) will vary about some value depending on the time 
of day, the season, etc. However, if the differences in z(t), u(t), 
and v(t) are taken between two different times, then it may turn out 
that the difference or change v(n) = v(t 4) - v(t) will have a zero 


mean, The general equation is then 


z2(n) = H u(n) + v(n) (3.16) 


wnere Ze = a(t) _ z(t) 
u(n) = u(t,,4) - alt.) 
v(n) = v(t 


oee - v(t) 





Agee 

It will turn out that there are other reasons for using the difference 
in variables rether than the actual values of the variables. The most 
significant is that using the actual values may result in subtracting 
variables whose magnitudes are very lerge compared to their difference. 

From now on, any reference to the previously defined symbols for 
voltage angles or power will refer to the changes of variable between 
two times rather than the actual values of the variables. Hence the tie 


line power flow model is 
PrL(n) = (YTL + YEOQX)DB(n) - PEQX(n) (35) 


and the boundary bus impedance model is 


DB(n) = ZaB(PEQT(n) + P8(n)| + ZBB PEQX(n) (3.18) 
where Pita) = PrL(t 4) - PIL(t ) 


rim) = fees) - PB(t ) 


DB(n) = D&(t,,,)'~ DaCt,) 


ie 


FEQI(n) = PEQI(t 4) = PEQI(t ) 


ene ae thee 


oO 
taj 
@>) 
rs 
—_ 
He 
~~ w 
| 


= PEOX(t .,) - FHOK(t ) 


for N sets of measurements n = 1,2,...,N. 
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eae TCR TONeOF LINEAR, STATIC SYSTEMS 


4.1 Conditions on the Svstem 

Both formulations for identifying the structure of the equivalent 
system have resulted in the same general equation, z(n) = H u(n) + v(n) 
where there are Nsets of input vectors, u(n), ane output vectors, 2(n). 
H is the structure matrix and v(n) is a pseudo noise vector. The 
general problem can be stated: 
Given: A set of input and output vectors u(n) and z(n), n =1,2,...,N 
Find: The structure matrix H and statistical information on v(n). 

Several assumptions on u(n), v(n), and H will te made and later 
these will be justified when the results are applied to the specific 
prodDlem. First elements of H are assumed to be independent. Specifically 
it is assumed that H is not symmetric. This is certainly not true for 
the specific problem, but various justifications will be given later. 
Second, v(n) is assumed to be a zero temporal mean process, or 
E v(n)t = Q. where Efat is the expected value of a. Third, v(n) is 
uncorrelated in time, or Efv(n) v'(m)} = 0 for n/m (n represents a 
sequence in time). Fourth, the covariance matrix for v(n) is assumed 
to be of the form TaN eh ONE E{y(n)x' (n)}} =n) = me) R. c(n) can 
be called a confidence coefficient for the ae measurement set and 
indicates the relative noise level for that set. Larger c(n)'s indicate 
less noise. Fifth, v(n) is uncorrelated with u(n), or Bfy(n)u! (n)t =e 

At this point there are many possible arproaches to the problem. 
Two will be used here. One is a tire tested mathematical approach 
which assumes a particular probability model for the noise, and the 


other is an engineering epproach, also time tested, wnich uses no model 





5s 
nor pememenemlity distzibution.of the noises It will turn out that 
both result in the same solution which will have a relatively simple 
form. The two approaches ere: 
1. Assume that the noise, v(n), has a Gaussian probability 
distribution (with a zero mean and covariance of ne Re 
Then find the maximum likelihood estimates of H and R. 
2. Assume no SG ced probability model for the noise and 
use the method of weighted least squares to fird f, Then 
make a reasonable estimate for R. 
The fact that the two approaches result in identical solutions is a 
well known consequence of using quecdratic criteria for optimization 
and estimation of linear systems. However, in acdition to analytical 
convenience, it is also well known that in many instances the Gaussian 
probability cistribution itself is a rather reasonable probebility 


model. 


4.2 Weximum Likelihood Identification 
First the solution will be found for the meximum likelihood approach 

assuming Gaussian noise. There are four known items: the set of input 

Peevomemii nine 1,2,...,N, the set of output vectors, z(n), n = 1,2,...,N, 

the corresponding confidence coefficients, c(n), n = 1,2,...,N, and the 

form of the probability distribution furction for v(n), n = 1,2,...,N, 

which is zero mean Gaussian, with a covariance matrix ey R. Note 

that R itself is not known. Therefore if estimates of H and R can be 


found, the problem is solved since knowledge of R completely specifies 


the statistics of v(n) when v(n) is Gaussian. 
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The probability distribution function for the Kxl vector v(n) is 


a a -1 
Spas la) Sa el 
p[yv(n)] = (Gow 7 *(n) RI] 2 e 2 cn n) 


The likelihood function of z(n) for a given H and R when 2(n) and u(n) 


are Kxl vectors and H and R are KxK matrices is 


1 
“2 oJ L2(n) .u(n) :H,R] 


p[2(n),u(n):#,R) = [(ar) ¢7*(n) RI] 


where 


i 


5[z(n),u(n) BR} = Sfe(n) - Buln} e(n) Rfz(n) - H u(n)] 


Since v(n) is uncorrelated in time (Efv(n)v' (m)} —JOefcrmna-unmtne 
joint probability distribution function for the set V = fv (1),v(2),...,¥(m} 


is the product of the distribution functions for v(1),v(2)...,v(Q). 


p[v(2)} p[v(2)] ... pfy(n)) 
1 


pv) 
N 


ON 5 e vw (n)e(n)R-+y(n) 
[(eny™ RI" ih o*(n)| x e ape, 


a 


NIH 


Iet p(Z:H,8) be the likelihood function of the sets 
Z = {2(2),2(2),...,2(m)}, and U = fu(1),u(2),... u(x} for a given H 


oma KR. Then 


1 
N > | 
P(e: Wen)” Wane Wc S(n)} e792 5:08) (4.1) 
n=l 
wnere 
N . 
3(2,U#R) = 2 Sen) - w wln)]' elm) R Yala) - # ula)] 


nu 


Tet A and R be the values of H and R which maximize the likelihood 
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function when Z and U sre the actuel measured sets. As developed in 


Apvencix D these maximum likelihood estimates are 


a N a 
H = 2 ¢ (jain )ut(n) i o(ndatoda (Lae 
n=) n=] 
and 
A 1 N A A t 
R = * > e(n)|z(n) -H u(n)| re -~H u(n)| (4.3) 


provided the vectors u(n), n =1,2,...,N are such that the inverse of 


N 
| 2 o(n)u(n)u' (n)| exists. 
m=) 

Two things should be noted. First, it is not necessary to know 


A 
R in order”to find This is significant in that it is not necessary 


A 

fe 
to solve for R and # simultsneously which would most likely result in 

A 
greatly increased computation requirements. The second is how R varies 
A A A 
with H. The quantity N]trace (2) is actually 2J(2,U:H,R) when R is set 
A 

equal to the identity matrix. &s seen in Appendix D, H is the value of 

H such that the matrix gradient of J with respect to H is zero. Then 


1 2 i jal : 
N34 [trace(a)] = $y era }=2 


A A 
or in some respect R is not sensitive to small errors in H. 





=i 
4.3 least Squares Identification 
The Gaussian distribution of v(n) was assumed in order to provide 
a probalistic framework, however the same equations can be obtained by 
assuming no probability model for v(n) and instead relying on the 
nearly universal metnod of weighted least squares. The sum of weighted 
squere errors between z(n) and H u(n) is 
we 
Be ym, \z(n) - 3 u(n)} e(n) T iz(n) - # u(n)} 
n=] 
where T is an unknown weighting matrix. Since S is the same as 
J(Z,U:H,R) when T = R71, from Appendix D it can be seen that the 
estimete for H did not depend on R anyway and the value of H which 
minimizes S for given sets 2(1),2(2),...,z(N) and u(1),u(2),...,u(N) 
will be the same as the maximum likelihood estimate for H. Under such 


circumstances, a reasonable estimate for E e(n)v(n)v! (n)} is 


N 
R= Efe(r)y(n)y' (nf my Z o(md¥ (nds (n) 


Ue 


i N A A t 
= e(e[ate) - Ba] [a(r) - 2 vl) 
ne 


| 
- 
€ 


which is the same es the maximum likelihood estimate equation (4.3). 





a Alpe 
4.4 Error Analysis 
The error analysis which follows only uses the covariance of v(n). 
For both the approach assuming Gaussian noise and the least squares 
approach, the covariance of v(n) was assumed to be E{y(n)v'(n)t = Ae Re 
Therefore results of the error analysis are applicable to both approaches. 


~ 


When the equation for z(n) is substituted into equation (4.2) 


nti n=) 


F N N N mi 
H = 2, c(n)H u(n)ut(n) + > o(abatobat (=) |] c(m)a(n)a (a) 
n=] n=] ne 
7 N N = 
H-Hs |: atoreonat |] e(ndu(n)ar (2) (Cae) 


One of the comitions on the solution was that v(n) was not correlated 
with u(n). It is interesting to determine the effect when that cordi- 
tion is not met. Let W be a weighted average of u(n)u'(n), n= 1,2,...,N 


N -1 [N 
wo | ote) \z c(n)a(n)a (a) : (5) 


me 


~ 


If v(n) is correlated with u(n) so that the weighted average of 


maui), o = 
+ N 
kt = 2; oP ie. e(n)v(n)u! | * z u(n)at(n)t CIS) 


A _ 
Then Ho eee CD 


and es N increases H =e wa inoveeO TO Zero as would be che cease if 
v(n) were not correlated with u(n). Error is introduced in the form 
of a bies which is proportional to the correletion of v(n) with u(n) 
and inversely proportional to the "power" of u(n). 


From equation (4.2) it may be noted that only the he element of 





SAS = 
z(n) affects the estimate of the ie element of E. Similarly from 
equation (4.4) only the ee element of v(n) affects the a row of 
the error ty - HK, This reflects the fact that E wes assured to be, 
in gereral, not symmetric, and that the ‘“" element of z(n) is affected 


' 
only by the 7” row of H and the ae element of v(n). Iet h, be the 


os mow of H. 


Then the estimate of the fee row of H is 


— 
— 


. N N ei 
h, = b oo) (oar on |] e(rdutoda' (aS 


ama DNe error as 


aA ‘ N N zal 
earls 2 clade (odar(n) | ctndu(ndst(n) (4.9) 


: ooh oun e 
The error covariance tetween elements of the i row amd j row of H 


is 





AIAN 


N == ili N 
=— [2 e(n)u(n)u! (n)| [= e(n)u(ndv, (n)] 
N N -] 
x > e(n)v, (oar) | Iz e(n)u(n)a (2) i 
n=] =] 


‘ -lTN oN 
7 2 o(n)u(n)u (n) | 2: 2. c(nda(n)a? (mde (nd=Gr, (adv, (| 


Na a 


N -1 
x |= c(n)u(n)u! cn) 


n=! 





Since Efy(n)y' (m)} =O for n Am am B§y (n)x" (n)} = stn) R 


A A ; N a 
Fi - ath, = hy) (a, = ny} Ee Ra \2 o(na(n)a (2)| (4.10) 


n=1 


let c be the average of the c(n) 's. 


i N 
ee >, ¢(n) 
n=] 


Then in terms of W, the weighted average of u(n)u'(n) defined by equation 


(4.5) the error covariance is 


R,. 
£ ii yl (ania 


&£S might be exrected, on the average, the error covseriance of the esti- 
inate is inversely proportional to the number of measurement sets. In 

e 2 a 6 -] 
terms of signal to noise ratios, the quantity c oe is a measure of 
the power of the noise v(n), and IN is a measure of the power of the 


Signal u(n). 
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See Oia GMICAMICN OF POAER SYSPEMS 


The problem of identifying the equivalent edmittance matrix and 
estimating statistics of the equivalent bus power injection vector 
for the external system model has been pleced into a general form, 
and the solution to the general form has been given subject to certain 
conditions. The two models for identifying the equivalent system are 


The Tie Line Power Flow Model 
PIL(n) = (YTL + YEQX)DB(n) - PEQX(n) 
The Boundary Bus Impedance Model 
DB(n) = ZBB[PEQI(n) + PB(n)\ + Z8B PEQX(n) 


Ween MW = lige 5 0+ 6 y Me 
In both models veriables such as DB(n) are the changes which occur 


between times t and t_. DB (n) = DB(t,) = rB(t). The general 


a 


problem is 


z2(n) =H u(n) + v(n) eee eee ly 





ee 


Ef v(n)y! (n)} = 7 


and the solution 


" N “N -1 

H = \2 ctodadnu tel] 2 e(n)u(n)ut(n)| 
n=] ==). 

A ie N 

R= y 2 o(n)[a(n) - Hun) [2(n) - E u(n)] 


was based on the assumptions 





ane 


l. Elements of H are independent of each other. 
> ESy(n)$ = 0 
BR, E§v(n) ut (my =© for nfm 





he Efy(n) vt(n)b = oy R= 1425-005 


5. Efy(n) ut(n)¢ 


iI 
10 


Now the effects of these essumptions on each model will be discussed. 
For both models if symmetry of H were taken into account at the 


beginning it would be necesssry to perform matrix inversions of dimen- 


1 


5 K(K+1) where K is the dimension of H. In developing the error 


Sion 

covariance equation for H it was noted that with H rot symmetric, each 

row of H is actually identified independently of the other rows.of H. 

With symmetry this is no longer true and the result is that all inde- 

pendent elements of EH (either upper or lower diagonal must be identified 

together. The size of the matrix inversion necessary can be significant 

as the dimension of H increases. For example, if the dimension of H 

is 10, the dimension of the matrix inversion is 10 when H is not 

assumed to be symmetric because es was shown earlier, each row of E is 

identified independently of the other rows and only one matrix is 

inverted. However, an inversion of dimension 55 is required when 

symmetry is accounted for since the 55 independent elements must be 

identified together. One way to handle symmetry at the enc is to 

treat the aa element and the fio element of tt es two estimates of 

the ce element of H (for i # j). Then using the corresponding error 
th 


variances, the weighted estimates can be combined. let a be the ij 


A 
element of ft and r its error variance, and b the ae element of H 
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and q its error variance, then the combined estimate c is 


ry [om 
+. 
QOlrY |} 2 lo 


il 
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Symmetry of ZBB and YEOX is merely the result of reciprocity since 
they represent passive networks (the sources, FEQT, P83, and PEQX have 
been moved outside). However, in the case of ZBB what is ultimately 
important is the product Z8B(FEQI + PB). By the principle of recip- 
rocity an input at terminal 1 causes an output at terminal 2 equal to the 
output at terminal 1 caused by the same input at terminal 2. But because 
of different demands throughout the network, the input at terminal 2 may 
never be as large as the input at terminal l. 

Fad the actual power injection and voltage angles been used, v(t) 
iar PEQX(t ) or ZBB FEQX (t_)) would most likely have a non zero mean 
value. However, because changes in values were used, it may be expected 
that v(n) will have a zero mean value provided al] the measurements are 
not taken when the entire system is moving in the same direction. For 
the case when the entire system is moving in the same direction, it may 
turn out that Bfy(n)y' (n-k)} =WOmecom < erceaucretnan | or 2, That ac, 
v(n) is only correlated with only the last one or two measurement sets. 
UiiesceoulasJeeaccountea for by a modification to the approach, but it 
is not anticipated that time correlation of v(n) will be a major source 
of error. Also, when the entire system is moving in one genera] direc- 
tion, Cepending on the sampling interval, there will still probably be 
some independent movement of certain buses in the opposite direction 


due to base loaded generators, load fluctuations, etc. 





- Ae 
The form of R(n) = E}y(n)v! (n)} Sa R was based on the 
= ei ene) - 
assumption that the covariance of the noise is dependent only on the 
time interval between measurement sets a(t), u(t.) and a(t i4)> u(t 4). 
It is then reasonable to assume that the covariance of v(n) = v(t .4)-v(t,) 


fs of the form R(n) = (t = t JR. There may be error in the meesure- 


ml] 
ments z(n) and u(n) but it was shown how such error could be lumped with 
Min). if v, (n) and v(m) are the errors of z(n) and u(n) respectively, 


then z2(n) = H u(n) + V(n) 


where v(n) = -v, (n) a n v fn) + v(n) 





Since v(n) is independent of v (n) and y(n), EYE (n)¥" (n}h = 9+ LR R 
' 

where Q = EAH vy, (7) = v, (n)] fy y(n) . v, (nyt and should be 

relatively constant since it depends on measurement error. The measure- 

ment error should be in the vicinity of one to ten percent while v(n) 

can be expected to be in the vicinity of 1C0% of H u(n) so the effect of 

Q on the estimate of R should be very minor. 

The last assumption, that v(n) is uncorrelated with u(n) can be a 
source of considerable error when this condition is not met. The error 
between fi and H was shown to be 

A N N -1 

H-H = =: etadutnda (| B e(n)u(n)at (n) | 

n= 


For the boundary bus impedance model 


DB(n) = ZBB]PEQI(n) + PB(n)} + 283 PEQX(n) 


v(n)u'(n) = 2BB PEQX(n)[PEQI(n) + PB(n)) 





= 2BB AQX|PX(n)PI'(n)AQI' + PK (n)FB'(n)] 


while for the tie line power flow model 





aoe 


Palin) =s2TL + YEOX)DB(n) - PEQX(n) 


a 


v(n)u'(n) = -PEOX(n)DB! (n) 


f 
-AQX PX(n)[AQI PI(n) + PB(n) + AOX PX(n)) ZeB 


= -AQX[PX(n)PI' (n)AQI' + PX(n)P3'(n)) zBB 


- BOX PX (n)PX' (n) AQX'ZBB 


where FQX and AOI were cefinec in equations (ry) anc’ (3-12). Both 
models have the term AQx] PX (n)PT! (n)AQt + PX (n) FB! (n)J in common. 
However, the tie line power flow model has the additional term 

AQX PX(n)PX'(n)4QX' = PEQX(n)PEOX'(n). Ordinarily it would be expected 
that a charge in power injection at one particuler bus an be indepen- 
dent of the changes in power injections at all other buses so that 
ESPX(n)PI'(n)} = Cee E{PX(n)F3"(n)} = 0. The exceptions might occur 
when the entire system is moving in the same direction as discussed 
earlier. However in the case of the tie line power flow model, 

E4 PEOK (n)PEQX' (n)} =R is not zero unless the exterral system bus 
injections cdo not change. This causes a bias in the estimate of H which 


depends on the magnitude of FEQX(n) in relation to FEQT(n) + FR(n), but 





in general it can be expected that the two will be of comparable magni- 
tude so that the bias introduced is in the vicinity of 100%. Using W 
and M as cefined by equations (4.5) and (4.6) 
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the bies in the estimate for the tie line power flow model is 


~H=MW wR ZB OW 


—— 
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hu> 


so that the estimation equations for a and R are really coupled for 
this model. 

The boundary bus impedance model and the tie line power flow model 
ere two specific approaches to the problem of identifying the equivalent 
external system model. In the former the inputs é@re Own System bus 
power injections and the outputs are bouncary bus voltage angles. In 
the latter the inputs sre boundary bus voltage angles while outputs are 
tie line power flows (which deperds on voltage angles). Eecause it is 
power which is Dought and sold, while voltage angles are cetermined by 
the line admittances and the distribution of bus power injections, it 
can be seen that bus power injections are the legitimate irndependent 
inputs and rot voltage angles. This is the cause of the bias in the 
estimate for the tie line power flow model. Both models need DB(n) but 
the tie line power flow model uses the tie line power flow vector, 
PIL(n), which in general will be more accurate a measurement than the 
estimate for the equivalent injections from IS, FEQI{n). It is possible 
to measure PIL(n) directly while the actual value of PEQI(n) is dependent 
on voltage magnituces and transmission line resistance eas shown by the 
nonlinear analysis in Appendix ©. However the solutions for f and R 
are more readily found for the boumary bus impedance model. 


Obtaining the estimate of YSOX from the estiméte of ZRB is straight- 





forwarc. However, if the estimate of YEQX is to ultimately be used in 
a linear load flow model for (CS, then as shown in Appendix B, the 


estimate for ZB3 can be used directly. If ZO is the bus impedance 





matrix for ©, 


DI zcs 
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6. SIFULATION RESULTS 


In order to test the method proposed, a computer simulation 
program was Ceveloped to generate data, identify the equivalent model, 
and check how well the identified model could predict line flow 
changes in 0S. Figure 6.1 gives a broad overview of how the simulation 
was conducted. First, using transmission line data, matrices to be 
used in the simulation were corputed. These incluee ZWS (impedance 


i 


Metrix for whole System), AO1, YIN ~~, (YEQI + YBOS), (YEQX + YTL), AOX, 


and YTL. For convenience, the confidence coefficients, c(1), c(2),..., 
c(N), were all unity. Matrices CHI and SGIN accumulated a running sum 
of z2(n)u'(n) and u(n)u'(n) respectively. In order to simulate the 
varistions in bus power injections throughout the network, P(n) (vector 
of bus vower injection chenges for WS) wes generated by reneom mmbers 


having an RMS (root mean square) value which was a fraction of the 


OQ 


nominal operating value of bus power injections Let One fraction, 
FCTCS, was used for 21] buses in © and another, PCTXS, was used for 

os e - @ “a mi 2 v 2 * * 
mie buses ineXS. (i.e. i PX’. (n) = \PCTKS x PX nom ) This main- 

i 

tained a relative scale for changes while at the same time it simulated 
independent changes in power injections. fs shown in Section 3, when 
the linear model is used, the bus voltege angle vector for % can be 


found using either 


Bien) PI(n) 
DB(n)| = ZiS | PB(n) (6.1) 


or equivalently 
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Figure 6.1 Basic Flow Diagram For Simulation 
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Ordinerily equation (6.2) or its nonlinear analog would be used to 
estimate the state of OS using real measurements of PI({n), P3(n), and 
PTL(n) (and any other redundant measurements). FEowever, for purposes 
of simulstion equetion (6.1) would have to be solved first to find 
DX(n) in order to generate PIL(n). Since either method results in 
the seme velues of DB(n) (in the linear model), equation (6.1) was 
used to simulete the output of a state estimation for DB(n). The 
simulated state estimation for OS results in errorless data. However 
as discussed earlier the only effect of measurement noise and state 
estimation soe is to modify the disturbence vector v(n). 

Everything described up to now was for the purpose of generating 
data for the identification routine. Using generated data, PI(n), 
PB(n), and DB(n), the input and output measurement vectors were formed, 
u(n) = AQI PI(n) + PB(n) and z(n) = DB(n). These values were also 
stored for later use in finding R. Running sums of z(n)u'(n) and 
u(n)u'(n) were maintained by CHI and SGIN so that at any time the 


1 eee 
. Using H and the 


current estimate H could be found by H = CHI SGINV™ 
stored meesurerent sets, en estimate B was fourd. Parallel to all of 
this, the actus] disturbance vector, v(n), found by v(n) = 2BB PEQX(n) 
= ZBX PX(n), was used to maintain a running sum of v(n)v'(n) in R. 
The estimated values i and R could then be compared to the actual 


values H = Z8BB and R. In orcer to compare how well the identified 


model could predict changes in line power flows, the estimate for 24%, 





ee 
the impedance matrix relating OS bus power changes to OS bus voltage 
angle changes was found by 


i aN 


y YI 0 AQI'| ZBB [AQr I) 
coe + (6.3) 
0 0 I 


fgain rancom variations of bus power injections were generated using a 


nominal base value and estimated bus voltage angle changes found by 


BY (n) PI(n) 
A | 
= 20S (Gay) 


63 (n) F3(n) 


were used to estimate line flow changes. The actual bus voltage angle 
changes using the same bus power changes were found from B(n) = ZS P(n) 
and were used to find actual line flow changes which could be compered 
to the values found using the identified model. 

Throughout the program, in order to present cata in a condensed 
form, for many veriables (z(n), u(n), v(n), actual lire flows, errors 
in vredicted line flows, and generated bus power changes) only the RM 
value, the mean, the minimum, and the maximum valves were maintained. 
It should also be pointed out that the basic approach is simple. Only 
a smell portion of the progrsm was involved in the actual identification 
of i and R. The rest is support for generating cate and evaluating 
results. 


A 
Using equation (4.11), the estimeted error covariance matrix, Paw 


@eeeoui of fois for c(n) = 1.0 





=5 55 
where 


1 1 N -1 | 
ane = i >, u(n)u! (n)| = N\ x Some = Noe SIGMA 
n=] 


Eence the estimated standard deviation of element ij of H is 

, x SIGNA,,)?. Since Ras should be somewhat proportional to the 
square of POTXS while aa should be inversely proportional to the 
square of POTCS, in the case of the lineer model with perfect measure- 
ments of z(n) and u(n), identificetion error is depencéent on the retio 
of PCTXS to PCTCS (percentage variation of XS bus powers to percentage 
variation of OS bus powers). 

One system presented in a paper [2 Jwas used here to test the 
method. The system shown in Figure (6.2) consists of 18 buses. For 
purposes of simulation Whole System consists of these 18 buses of which 
buses 1,2,3,4, and 18 belong to IS, 5 through 7 belong to BS, and 8 
through 17 belong to XS. Bus 18 is the reference bus. Other parti- 
tionings could be used but this one does result in matrices (YEQX + YTL) 
ana (YEQI + YB°S) which sre comparable in size. This avoids a case in 


which coupling of BS with IS is very much stronger than BS with XS. 


For this system (YEQI + Y8CS) and (YEQX + YTL) are 


3463 aeO.0 ~4.07 
YEOT + YBOS = 0.0 OTe eS 
=7,07 ~=9.42 | 65.9 


Mies -Ont -5.36 
Mie + YL = | -6.1% 14.2 -8,12 
aoe one he 12.5 
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Both the effects of increasing disturbances anc increasing 
measurerents are sirulated here. In Figure 6.2 and Table 6.1 are the 
results of increesing disturbances (PX(n)). The RMS of bus power 
changes for OS in each case was 10% of nominal while RMS valves of 
lemon) were O07, 1.25%, 2.57, 5.0%, 10%, and 2C% of ee Another way 
teesey this is that "signal Mepnorse revicos were ctl, G:1, 4:1, 2:1, 
l:1, end 1:2 respectively. The same sequence of rancom numbers were 
generated each time. Included in Table 6.1 ere actual values of ZBR, 
actual errors in identifying it, and the estimated standerd ceviations 

A 


using R in equation (4.11). Also included are the values of 


“N 
(YEQX + YIL) estimated from ZBB, and the errors of predicting changes in 





line flows using 2B. In a similar format, the results of increasing 
the number of measurements are summarized in Figure 6.4 and Table 6.2. 
Changes in both OS and XS bus vower injectiors had RMS values of 10% 
of the nominal operating voint while a was identified using 16, 32, 64, 
128, 256, and 512 measurements. 

The results show eee rat only is identification error commersurate 
with the noise and number of reesurements, but also estimeted standard 
deviation gives an accurete measure of the identification error. 

Because the same random numbers were generated for each value of PCTXS 
in the first cease, the actual errors shown in Figure 6.2 are nearly 
exactly linear since as was shown earlier, for equal numcers of 
measurements identification error Cepends on the ratio of PCTXS to 
PCTOS. For the case of increasing measurements in Figure 6.4, ZBB 


Jel 


does not have a steadily cecreasing error es with ZBB.., but it is still 


22 


in the range of the estimated standard ceviation. 
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To contrast identification mofels, in Table 6.23 the results of 


aN 
identifying YEOX + YIL by the tie line power flow model along with the 


a 
estimsetec stancard deviations are compared with YEQ* + YI as found 


centified by the boundary ous impedarce model 


pe 


nN 
from 223 which was 
using the sare cata. RBesices idcentifving the wrong values, the tie 
line power flow model gives inaccurate estimetes of error. 
In Figure 6.5 is another network presented in a paper by Stage see 
%S consists of the higher voltege network and ¥5S is the lower voltage 
network. The "tie lires" ere really transforrers. Using the boundary 
bus impedance model an equivalent network is identified to replace the 
21 buses of XS, so that can be analyzed with 9 buses instead of the 
original 30. The results in Table 6.4 again show close correspondence 
between ectual error and estimated standard deviation. 
The simulation carried out was certainly not fully realistic. 
memever the realistic aspects include 
1. Values of bus power changes were found first, then a state 
estimation for OS was simulated and sopropriate values were 
pessed to the identification portion of the progrem. 
2. The effects of increasing meesurements were studied under 
the assumption that XS would vary es much as OS. 
The simulation was unrealistic in that 
1. The nonlirear model was not used to generate data with the 
result that a linear identification model was used to find 
paremeters to a linear system. To gererate data from a 
nonlirear model properly would require a complete load flow 
for each measurement set. (It would be sirple to vary voltage 


magnitude and voltage angle end then find bus power, but bus 
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power must be the specified variable as it is the incepencent 
input, not voltege angle.) 

2. Bus power injections were indepencently generated whereas 
realistically they will probably be correlated in space end 
time. There were no comitions on the solution for correla- 
tion in spece. In fact FEQX(n) itself hes a covariance matrix 
E $ FEQX (n)FEOX' (n) 5 which is in general not diagonal. However, 
the solution did assume no correlation in time. fs discussed 
earlier in Section 5, this effect can be accounted for by 
modifving the epproech. 

Based on the results of these and other simulations certain con- 
clusions can be made about the boundary bus impedance model. First, 
with enough meésurements it is possible to identify the equivalent 
mocel under conditions of the simuletion (independent movement of the 
bus injections). The number of reasurements may seem rather large but 
of course the noise is @lso rather lerge. The results show that as 
expectec either a decreese in the cisturbances by a half or using four 
times es many measurements éepproximately halves the error. Secord and 
vossibly even more important, the error vredicted by using the estimated 
R in equation (4.11) egrees very well with the actual answer. This 
means that an accurate estirete can be predicted from only the measure- 

A A 
ments used to ifentify H so thet the limits to which HEH can be trusted are 
also known. This is of course not the case with the tie line power flow 
model. 

A very irportant question which may arise is how much data would 
be required to identify the system. The simulation was carried out 


under idealized circumstances and as was shown earlier, for the 





wea 
Simulation the icentification error for equal numbers of measurements 
Gepencs only on the ratio of XS bus power changes to OS bus power 
changes. In ectuality however there are upper end lower limits. If 
data were used from samples taken too close together in time so that 
changes were smell, the effect of measurement noise and state estima- 
tion error could become very significant. On the other hand, the 
System can vary over only much of a range, so thet if the time 
between measurements were too large, then because of cyclical patterns 
of power requirements, there might not be much new informétion after 
the first few measurements. To determine the emount of data required 
anc the best interval between samples would require a study using 


actual operating cata. 
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YEOX “ew 
Actual Value 1S -6.14 -5.36 
Error BBIM ECO -1, 27 1,30 
Error TLPFM - 20 = a, 10.66 
Reo. 1 Lem 519 2.0 39 
Actual Value -6.14 17726 -8,12 
irror EBIM 07 -1.87 2.97 
Error TLPFM Le52 -19, 37 14.18 
Bao... ligeM we & oe ol 
Aetual Value ~5. 36 -8,12 12.47 
Error 28™ .52 ee 2.79 
Error TLPFM Saal -22,24 6. o 7 
Mo... Lee 60 1.29 1,26 


Actual Value -- Actual Value of YEQX + YTL 


7 
Error BBY -- Error in Identifying YEQX + YIL using the Boundary 
“™, “N 
Bus Impedance Model YEQX + YTL = ZBB - (YEQI + YROS) 


vo 
Error TLPFM -- Error in Icentifying YEQX + YUL using the Tie Line 


Power Flow Model 


A 
E.S.0. TLPFM -- Estimated Standarc Peviation using R from the Tie Line 
Power Flow Model 


OS and XS varied equal emounts (1C% of nominal) and the same 


256 measurements used for both models. 


Table 6.3 Comparison of Identification Models 





oe 


285 for 20 Bus System 


Ectual Value BGGAY, C459 
Eetual Error -.CC16 .CG02 
mes Ds blo C020 
Eetual Value .OL89 C665 
Actued Brror -.COl9 2@0C7 
Seb, 0621 0022 
Actual Velue C291 .C663 
Actual Error -. C02] C010 
See C025 0026 
Bw... -- Estimated Standard Deviation using 2, 
7 > 
Z38 Identified by the Boundary Bus Impedance Model 
1Cé variation of bus povers in both XS and CS. 


256 measurement sets used. 


Table 6.4 Identification of a Lower Voltage 
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C179 
sO2 le, 


.C663 
.0O184 
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7, USE IN POWPR SYSTEMS--SOME PRACTICAL 45PECTS 


I~ 
} 


: Use of Bus Fower Jnputs 
From equation (4.10) which relates the error covariarce of the 

estimate, it is obvious that the input directly affects the error. 
For the boundary bus impedance model u(n) = PEQI(n) + PB(n). As 
discussed ctreviously and shown in Appendix C, FEQI(n) really depends 
on the operating point in a nonlinear manner so that it is possible for 
FEQI(n) to introduce significant error. The error from P8(n) on the 
other hand is only due to measurement noise. Yor this reason, it 
should turn out that the identification method presented will be most 
successful when the bourmary buses have relatively large power injec- 
tions (generators or loads) which supply much of the known variation. 
Oe other words, it would be nice to have good, strong input signals. 
In actuality IS may be of reletively large size so that contributions 

f PI(n) to PegI(n) = AQT 


T(n) from far ecross the network may be of 





questionable value and accuracy Cue to the linear approximation of AQT. 
Since intuitively, buses in I closest to BS will have the most affect 


on PEQI(n) and involve the leest error, it may be prectical to use only 





part of the contribution of PI(n) to PEQI(n) ard lump the rest with 


PEOX (n). 


1.2 .ModeleVenification 
Although contrary arguments could be mace, the boundary bus 


impedance mode] is probably the better of the two models since it does 


not involve estimating a significant bias. Fowever, the tie line power 
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flow model is still valid for verification purposes. For the bourdary 


bus impedance model from equation (4.3) 





A 1 N A A 
| > ¢(n)¥(n)¥! (n) 
n=l 
N wea ax 
. Z, e(n)ZB8 PEQX (n)PEOX! (n)ZBB a 
nS 


il 


A as 
If R is premultiplied anc post multiplied by ZBB ~ then an estimate of - 


the covariance of FEOX(n) is 

A A LZ Nwcta 

R mae S255) 3 oh (72) 
Now if measurements of PTL(n) were taken at the seme times as measure-~ 


ments of DB(n) and PEQI(n) + FB(n), and if the estimate of YEQX is 





1 


ao “Ni 
YEQX = ZBR™* - YEQI - YBB 


Zax A 
then equation (7.3) evaluated for YEQX should be close to Reeox of. 


eqmetion (7.2). 


N 
>, ¢(n)PEQX (n)PEOX ' (n) 
ned 


ro | 


1 N 
-is { [et (n) - (YPL + YEQX)BB(n)} 
nel 


x [PDL(m) - (WEL + yBox)Da(n)}" e(n)t (7.3) 


If the two were not close, then the tie line power flow model would be 
indicating a leck of validity for the estimate of Z23B from the boundary 
bus impedance model. The tie line power flow model could also be used 
"on-line" once YEQX hed been identified to provide a continuing check of 
the validity of the identified equivalent system. Should the measured 


ao 
value of PTL(n) vary from (YITL + YEQX)DB(n) an amount not commensurate 





ey 


; ) A ee: : 
with ae ReEox? that would indicate a change in XS had taken plece. 
The change of course could either be in transmission line status 


affecting YEQX or a significant change in PX(n) affecting PEQX(n). 





7.3 Estimating Equivalent Power Injections 

Once veox has been identified, the actual values of FEQX (t, ) could 
be estimated. With enough measurements it might be possible to develop 
an approximate daily pattern for PEOQX(t) so that they could be used as 
pseudo measurements [8] for the purpose of predicting a future power 
flow situation. If such were the case, the equivalent system could be 
extended in use from predicting changes in @ curing which the change 
PEQX(n) would be small to use in predictirg actual line power flows for 


given power injections in G. 


7.4 Inputs Which Decrease Identification Error 

It has been assumed that the identification would be primarily 
passive in nature (i.e. by "listening" to the system) so that any power 
changes or voltage ancle changes used to identify the equivalent system 
would be due to changes in the system's load demands. It was also 
argued that se R was a reasonable approximation to the noise 
covariance where c(n) is inversely proportional to the time, | - Oa 
between measurements. However, any large, relatively fast changes in 
OS would provide an input, output measurement set with relatively little 


noise (PEQX(n)). These changes could be planned or accidental. If 


planned changes were used, it would be desirable to obtain a maximum of 





Te 


invormation with a minimum of change. For this case it might be 
possible to formulate the problem in terms of minimizing a cost 
function subject to a constraint such as finding the inputs 


u(n) (n=1,2,...,N) which minimizes the cost 
A 
c = f, (error in H) + f,{u(n)u'(n), (ae = Ieee | 


subject to 


z(n) = H u(n) + v(n) (n=. 2eemes ND 


where fs and fy sre scaler functions. 


7.5 Use of the Faquivalent System When the Fxternal System is Known 
Aithough the linear load flow model was used to solve the problem 
of finding equivalent systems with no prior knowledge of the external 
systers, there is still an advantage to using the model when enough is 
known about the external system so that, for instance YEQX and PEOX 
could be found analytically for a particular system. The linear model 
clears the smoke so to speak, so that a good approximate victure can 
be easily obtained. For exemple AQX (and AQT) gives a general overview 


of how power is routed by the network while a comparison of YEQX with 





YBB shows at a glance the relative coupling between the two networks. 





yy 

7.6 Identification of the Nonlinear Model 

Use of the linear model may be sufficient to predict line power 
flows amd voltage angles. However, in order to predict the effect of 
XS on voltage magnitudes within OS, the full ronline@r model is neces- 
sary. Peveloped in Appendix C are the equations releting complex bus 
voltages to complex bus power and current injections. The computation 
required anc complexity of identifying the nonlinear equivalent model 
would likely be very much more than for the linesr model. The increased 
accurecy in predicting real power flows would most likely be minimal, 
but the linear model cannot predict i pcan Gue to the struc- 


twee of XS. 
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8. CONCLUS IOts 


&,1 Areas for Further Study 

The tie line power flow model produces a biesed estimate of the 
equivalent model parameters, however, if an unbiased estimate could be 
obtained so that accuracy was comparable to the boundary bus impedance 
model, there would probably be an overwhelming advantage to using the 
tie line power flow model. One of the most important reasons is that 
it uses variables from a local region only. PIL(t) is directly 
measurable and in order to find DB(t) state estimation would be required 
only in a small area around Boundary System. (Power coming in from the 
rest of Own System could be treated the same as PIL(t) for purposes of 
state estimation.) Also the linear tie line power flow model appears 
to be less of an approximation to its nonlinear analog (see Appendix C) 
than the linear boundery bus impedance model is to its nonlinear analog. 
Therefore, it is suggested that a study to find a method for eliminating 
bias in the tie line power flow model would be worthwhile. 

The simulation carried out used a linear load flow model for the 
purvose of generating data used to identify the equivalent model. This 
Simplified the simulation and it also served to eliminate errors which 
would be caused by variations from ideal conditions (low line resistane 
and voltage magnitude nearly constant). A more realistic simulation 
should be made to find how cata from a nonlinear model affects accuracy 
of the identified linear model parameters. Also the net power flow summed 
over all the tie lines is usually held to a scheduled value. The effect 
ef this on the identification should also be studied. 


As mentioned previously, the identification accuracy could be 
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increased by varying the generation pattern for the purpose of probing 
the external system. However it would be necessary to maximize informa- 
tion and minimize the variation required. It is felt that one step in 
the proper direction might be to move one generator at a time since this 
would eliminate eny interference among the inputs. A look at the boun- 
dary bus impecarce model written in terms of the general model 

z(n) = H u(n) + v(n) shows that if elements of u(n) were all zero except 
for element i, then the io) conkemayort E would be en {2 (n) =e Cay 

A study into this problem should result in a more a eee solution 
to a rather complex probdlem. 

An obvious extension of identification of the linear model is the 
identification of the full nonlinear model. It is felt that the con- 
plexities involved would be an order of magnitude above that required 
for identification of the linear model. The accuracy gained in predicting 
voltage phase angle changes would probably be small, but voltage magnitude 
changes which have been ignored, would be available from a nonlinear 
model. If a search type of solution were usei, then of course the 


solution to the linear model identification would give a starting point. 


8.2 Summary 


The problem of predicting how an electric power network will interact 
with other networks has been solved by reducing the entire outsice world 
down to a relatively small equivalent model. The problem is complicated 
by the facts that (1.) no knowledge of parameters or variables of the 
outside world can be assumed, (2.) identification must be primarily 


passive, and (3.) "signal to noise ratios" can be expected to be about 1:1. 


~20- 


& linear model for the overall system was developed and its use 
justified with respect to the more conventional nonlinear model by con- 
sidering accuracy versus computation requirements, insight, and com- 
plexity of the identification approach required. Two models were then 
proposed with each having certain advantages amd shortcomings. The tie 
line power flow model uses data which is relstively accurate and locally 
available, but it produces a biased estimate. On the other hand the 
bouncery bus impedance model produces an unbiased estimate but it 
requires cata from all over the system and also requires manipulation 
of the equations for own system. It was argued that changes which occur 
in power enc bus voltage angle between two times should be used with the 
hope of settee zero mean disturbance and also to avoid the use of 
matrices wnose elements have magnitudes very large compared to the 
determinant of the matrix. Since both models were in the general form 
of a linear, static, input-output relation with additive Als tucbenecee 
a solution to the general problem was found under assumed conditions. 
Then the solution was applied to the two models and it was shown how 
the tie lire power flow model produces a biased estimate. The simulation 
made to test the method was ciscussed showing the excellent correspon- 
dence between theory ard simulation results. Finally, practical aspects 
of the problem were briefly discussed including model verification, 
estimation of the cisturbances, system probing, and use of the linear 
model reduction method for purposes other than identification. 

This method is proposed as a means of modeling the outside world 
and fincing the parameters of that mcdel.. The actual solution is ina 


very simple form and could probably be implemerted as an extension of 
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state estimation techniques. It is by no means the ideal answer to how 
the outside world should be modeled, but it does allow a prediction of 


how the outside world will interact with the system of interest. 
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APPR IX A. NIWLIRBER LOED FLOY EQUATIONS 


For the pi transmission lire model shown in Figure 2.1], the current 


flowing into the line from bus i is 


ie.) = E. yS., + (<8 


ik 17 Ey) Vay (A.1) 


when the line is connected tetween buses i and k. aap is the shunt 


capacitive admittance, Ss is the series admittance, ard a and Ey are 


the voltages at buses i and k respectively. A411 are complex quantities 


and 
é 
aoa faves ee 
ie ec a JO a | (A.2) 
_ ~jS. 
a deka 1 
‘ee 7 2 + (4.3) 
ik ik 
= JD; | 
ee ; (A.4) 
a JD, 
EY = EY e (A.5) 


where R LY and yS.) are the resistance, reactance, and capacitive 


DE acs k 


susceptance of the transmission line between buses i and k. D; and Dy. 
are the voltage ohase angles at buses i and k as measured with respect 


to a reference. Equation (A.1) can be rewritten 


SEY Saga V5 >) 
By Kirchhoff's current law, the sum of currents entering the transmission 
lines conrected to bus i must equal the current being injected into the 
bus by gererators or loads (positive for generators and negative for 


loads). uy is the current injected into bus i. 
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as 2. a Z WR t Fy) > 2 ys (4.7) 
k k k 

Let ee be a complex bus admittance matrix whose diagonal element in 

row iis the sim of shunt and series acdmittances of 211 lines connected 


to bus i 


> @ Gaae ae ee (A.8) 
k 


and whose off diegonal ik and ki elements are the regative value of the 
series admittance of the line between buses i ard k (or zero if there is 


no line between buses i and k). 


CTene “aie ee Tie | (4.9) 


i is a symmetric matrix. Then for complex vectors is and 


=bu's ae 


whose eh elements are the current injected into bus i and the voltage 


at bus i respectively, equation (A.7) can be written in matrix form as 


oe : ae ae (A. 10) 


Wnen Kirchhoff's current lew is epplied to the entire network, it is 
obvious that for a system with K+l1 buses, the current injectiors into 
K of the buses determine the current injected at the (Kei) bus. 
Similarly, the voltage angles sre relative to one another so ore bus 
must be designated the reference bus for voltage angle whereas the 
reference for voltage magnitude is ground. 

Real and reactive power injected into bus i, Pp. t JQ; in terms ot 


the current injected and bus voltage is 
P, + jQ, =E, I (A.11) 


or Deedee = E (A.12) 
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where 3a is the complex conjugate of voltage at bus i. Using 


equation (4.7) along with (4.8) and (4.9) where Y., is defined, 


ik 
equation (4.12) is 
sae 5 7 
Fo. = jO. = 3: c &. 
st a i Kk @ ak 
-jD, Kl 4D -je. 
= ee i Ey e : Yep e - (A.13) 
i a i 
os = 1k : . Be 

ces a ele (2.14) 


The load flow equation, (A.14), can be used for networks consisting of 
transmission lines end fixed tap trensformers. However for other cir- 
cumstances such es voltage controlled buses, tap changing trensformers, 
and phase shifting transformers the equation must be modified. Stagg 


ard El-ébiac gle) Fe one reference for these anda other variations. 
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NWPPFEIDIX 8B. LEMEAR LOAD FLOW EQUATIONS 


The linesr approximation for real power flow in a transmission 


line between bus i and bus k is 


Cu 


Pi, - Bee. Om DD (Cr) 
where Dp. and De are the voltage phase angles at buses i and k as 
measured with respect to a reference, E, and FE, are the nominal values | 


i k 
of the bus voltage magnitudes, and Yay is the transmission line edmit- 
tance magnituce. This approximation is subject to the conditions 
lt; De = Beistsrais. 
at K 
2. The bus voltage megnitudes do not vary much from their 
norinal values. 
32. The ratio of transmission line resistance to reactance is 
Sali, 


- 


ww 
For convenience it will be assumed that E = EL = 1.0 (or equivalently 


Vax is normalized for i and E) so that equation i) dis 


Phy = Vy (Py - Py) = Pavan - YayPy (3.2) 


Assuming no losses within the bus, by continuity of power the power 
being injected into a bus (positive for generators and negative for 
losds) is equal to the algebraic sum of the power flows into the trans- 
mission lines conrected to that bus. Let P, be the power injection at 


bus i, then 


a a a Z Yan 7 2 ener (B.3) 


Define P to be a vector whose oe element is the power injected into bus 


i, and define D to be a vector of corresponding bus voltage angles. 





Then relation (B.3) becomes 


y= (ea) 


By comparison with (B.3) it can be seen thet the elements of Y are such 
that the diagonal element of row i is the sum of acmittance magnitudes 
of the transmission lines which are connected to bus i, Ye; = >, Vix 
k 
The off diagonal elements ik and ki are negative value of the acmittance 
eG } 3 } } ° — x a wee e 
megnitudce of the line between bus i and bus k, Tey 4 Vsx 


Therefore the matrix Y is of the form 


% Vix “Yy2 “Y73 ~S= 
“Yy0 ZY, 493 + 

Y = (B.5) 
gals cee Zs ¥3 i 


Because each bus is usually only conrected to a few of all the other 
buses, Y is normally a sparce matrix. It is readily observable that 

the sum of all elements of any row or column of Y is zero. The bus 
voltage phase angles must be measured with respect to a common reference 
so ore bus must be designated as the reference bus where voltage angle 
is screcified. Also, because the linear model is lossless, by continuity 
of power, the algebraic sum of elements of P is zero. Hence one row of 
equation (8.4) is redundant. If bus i is chosen as reference for vol- 
tage phase angle, then let ie be the vector D with element i deleted, 


Het °F . be the vector P with element i deleted, and let } 


aw 


matrix Y with row i and column i deleted. Then equation (B.4) is 


has = be the 


P on D (B.6) 


=OUs =—pus =bus 





OT. 


rus ang bus 


System (IS), FB and DB for Roundary System (BS), and PX and DX for 


can be separated into vectors PI and DI for Internal 


External System (XS). Similarly separating Y into submatrices 


bus 
PI Gory oe |) as 
pe) | eee ae aes) ash ||| sas (3.7) 
PX Mr vee VX¥X DX 


However by Gefinition there are no transmission lines between buses of 


Peranc AS so VIX = 0, 


PI Wis YB a | for 
ra, = | ym yep = yax| | pa | (B.8) 


PX 0 yBX' YX | | DX 


The inverse of the bus admittance matrix Y is the bus impedance 


Mc 
matrix Zs such that 
ae : ae ere (B.9) 
or DI 211 ZB 20a BX 
Dee ie 283 Z2K PE (2 alg) 
Z in terms of submatrices of Y can be found by solving 
—bus —bus = 
Z,. Y = J for 2 Voene Lic) tie aGent1u matrix. 
molice ~OUus = = aS ra ‘ 
oe 6 2 a BO Io Q 
Zig" ZBB8 ZBX| |YB' YBB YBX| =/0 JI 0 (3.11) 
Ze ZBX ! 2XX 0 Moe YOR OO of 


—e= =~ — 





~se~ 
Multiplying the middle row of ee by the left column of ae (nere 


rows em columns will refer to row and column submatrices.) 


: 
< 


IT+ Z3B YIB' = Q (eae 


—— = ——— 


ame solving form ZIB' 





Z1B' = -288 YIB' YIT* (3.13) 
Multiplying the middle row of ae OVA eer I eneecos aan of Lee and 
solving for ZBX 

223 YEX + ZBX YXK = O eas} 

ZBX = -Z3B YRX Yxx"* (B.15) 
Multivlying the middle row of 2 by the middle column of Y and 

; bus bus 
solving for 288 with the use of (B.13) and (8.15) 
aap" VIB + 208 Wee + 2k YPA' = J (Beee) 
-233 YIB' YIT* yIB + 233 YBB - 283 ¥BX yxx7yBx" = I (B.17) 
-1. . cai -1 
ZEB = (ee Ye ek EX EK") (B.18) 


Multiplying the top row of ae byevse lene column-o! Y¥ and solving 


—bdus 
for ZII with the use of (2.13) 


ZIT YII+ ZB YB! =] (8.19) 
7 ee ese iT (3.20) 
gtr = vit? + vit yp ges ypBt yir? (B.21) 


Multiplying the top row of ae by the right column of yy and solving 


us 


for ZIX with the use of (B.13) again 


Z1B YBX + 21K YXX = Q oe; 





fe 


ER 


= ~ 213 vax yxx7t 


ee 


ry, Tee VEX 


ea mee eee ee ee 


it 


| 
I< 


Tet matrices YEQI, YEQX, 491, and AQK be defined by 


Q 


ta 


_ 


| 


YEQX 


iS 


o> 
© 
Ps 


Then equations 


If i i B 


ie. 
io) 


il 


i 


= ae? as 


al 


-YBX YXX"~ YBX! 


= -YB' yu 
for ZEB, ZII, ZB, ZX, ard Z28X can be written 


= (YEQI + YBB + YEQX)™ 


eee ee eee 
= AQT! ZEB 


I 
> 
Oo 
ee 

ts 
WwW 

In 

2 
Ps 


it 

Ww 
WJ 
to 
> 
a) 
S 


Using (B.1C), the equations for DI and DB are 


Zi 


hal 
N 
od 


PI Z| PX 


—_ 
= CT 2s » ance ——=- > = 


or using (3.29) through (B.30) 


(Beo3} 


(B.24) 


(B.25) 
(B.26) 
(B.27) 


(B. 28) 


(B.29) 
(B.30) 
(B.31) 


(Bae2) 


-_— 
Ww 
A) 
8) 
ae 


(B.24) 





DI err 
DE 0 
a 
i 


we 


0 AQI'| zee [sagt I}) [PI 
+ 

0 I PB 

ZBE AQX PX 


| 


From equation (B.35) the key role of ZBB is csvious. 


(B.35) 
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PPEBWDIX C. TOWEIMIAR MODEL IDM Mr ICAT ION 


To avoid notation clutter in this appendix, all voltage and current 
vectors and acmittance matrices will be understood to be complex vectors 
and matrices even though there is no bar over the symbol (for a complex 
quantity) or bar under the symbol (for a vector or matrix). For example 
E and I are complex vectors. 

The complex matrix equation for relatirg bus voltages to bus cur- 


rent injections is 


ae = ee Eis (C1) 


— and Ens can be Civided into separate vectors for Internal Systen 
(II and EI respectively), for Boundary System (IB and EB), and for 


External System (IX and EX). Dividing Y into the corresponding 


bus 


submatrices and noting that by definition of IS, BS, and XS there is no 


immediate coupling between IS and XS, equation (C.1) can be written 


ie mor =~ «VIB 0 EI 
ml -=|yYm' yYeB  yBx| | EB (C.2) 
x O voc eee) wes 


The only coupling between OS ard XS is due to YBB which can be separated 
into two matrices. That belonging to CS will te designated YBOS, and 
that belonging to XS will be designated YIL, so that YEB = YBOS + YTL. 


Separating equation (C.2) 





Soo 





6 EI 0 0 0 EI 
TB = Pet 4 1tee Va YRY. EB (C33) 
1 EX 0 vext “YXX ax 
Removing the 
II Water YIR EI 0 
= + (Gav 
TR vy’ ss yaas| | ER YTL EB + YBX EX 


The quentity YTL EB + YEX EX is a vector of line current leaving the 
buses for BS over the tie lines between BS and XS. Designating this 


the tie line vector IIfL 


Tiga ‘ar ie ll 
= (gS) 
T8-ITL YIB'  yscs EB 
TL = YTL EB + YBX =X (RS), 
Removing the bottom row from equation (aE 
IX = YBX' EB + YXX EX (C. 79 
ana solving for BA 
ex = -yxx7) yex! eB + yxx7? x (C.8) 
Substituting equation (C.8) into (C.6) 
IL = YTLEB - yBX yxx7? yaxt EB + yex yxx7? x (c.9) 


Define an equivalent external bus admittance matrix YEQX and an equi- 


valent external complex bus current injection vector IEQK by 





Oss 


YROx = ~YBK YXX ~ YSBX! (Ge 10) 
TEQX = AQXY IX (Capa) 
baie LOX = -YBX yxx72 (C.12) 


Then equation (C.6) becomes 
ITL = (TL + YEQ@ap - GEex (C.13) 


and if equation (C.13) is substituted into equation (C.5), since 


Mone ="yo05 4° VTE 


pels Yoh ais EI 
> (C.14) 

IB8+ IEQT Yi! (YBB+YEQX ) EB 
In the same way, an equivalent complex bus admittance matrix am an 
equivalent complex bus current injection vector for IS as seen by BS 
can be found. The top row of equation (C.14) is 2 

II = YII EI + YB EB (Creep 
Solving for EI 

zs peer | -1 

Ele=-Yrr Vip ED 4 auioeeem (C236) 


and substituting into the bottom row of equation (C.14) 


a a 


18 + IEQI = -YIB' YII~ YS EB + YIB' YII- II + (YBB + YEQX)ES (Cig) 


Defining an equivalent internal complex bus admittance matrix YEQI, and 
an equivalent internal complex bus current injection vector JIEQI by 


5 


VEOL = ib! YIlo= vas (CsI) 
mofe-ere il TE (C.19) 

= 
where hoe eey Tp! YLT (C.20) 





equation (C.17) becomes 
TEQI + IB + JEQX = (YEQI + YBB + YEQX)EB (6.21) 


Just as YBB was divided into YTL and YBCS, so also YBCOS can be divided 
into that belonging to IS designated as YEBS, with Y3C0S = Y8IS + YBBS. 


Then equation (C.21) is 
T2QI + IB + TEQX = (YBIS + YEQI)EB + YBBS EB + (YTL + YEQX)EB (C.29) 
Haus far equations) (6.5 ),0 (6.12) Ge ee ome. 21) earrentiec 
Pertinent equations in terms of bus current injections. To place them 


in terms of bus power injections let a diegonal matrix whose elements 


are the complex conjugate values of bus voltages be 


x 
ES Co 
Xa 2 
[eo] 0 E, 0 -- 
0 Cael cee ; 


OS aw at aw & ae ow Fae 8D aw oe =] av 28 oe oF == aw 


(Note that [E) is a diagonal matrix whose diagonal elements are iden- 
tical to the elements of the vector E.) Also, let PI and QI, PB and 
Q8, and PX and QX be the real and reactive bus power injections of 


IS, BS, and XS, Then because for any bus i, P, - iQ, = E, I, 


ee cree ers (C225) 
PB - 308 = [FB | B (C.24) 
Poe on = ee x (Ces) 
pre - jqtu = (ss J mu (C.25a) 


or 





SoG, 


axel 
m= fer] (pr - jaz) (C.26) 


«cl 
X sain ] Gee Foo C2) 
2 

Multiplying equation (C.21) by [EB J 

% % 
[e8 |\(12QI1 + 1+ ax) = [EB )(YEQI + YBB + YEXX)EB (C28) 
Define equivalent internal and external bus power injection vectors by 

% v7 
eer \eraier leer =o) (Gea) 


PEQI - jQeQqr = [eB] rar 


il 
li 


x x x47 
PEQX - jgeqx = [es ] rox = [eB } aoxfex ] (px - jax) (C.30) 


Then (C.28) becomes 


(PEQI - jQEQI) + (FB - 308) + (PEQX - jqzQx) = [eB] (YEQI + YBB + YEQX)EB 
(C.31) 
For the following a bar over a symbol denotes a complex quantity 
while matrices and vectors are denoted by underlined symbols. Using 
equations (C.5), (C.22), (C.24), and (C.25a) the nonlirear load flow 


equation for (CS is 


PI - jQI fe ae iene ley |) || ea 
2 (Gag2) 

ace Ne amet pee <_< 

(2 = 968) - (BRE 9971) o {zB | ym yses| | BB 


—— 


which, if YEQX and PEQX - jQEQX are known, becomes 





PB - jQ8) + (PEOX - jQE0x) o [rE] 


al ESI 


(od | 


~96- 


‘OI ltt] o 





The nonlinear analog of the linear load flow tie line power flow model 


for identifving YEOQX is 


PIL - jorL = [EB } (YTL + YEOX)EB - (PEOX - jQEQX) 


el ae Cm ore oe = eee 


and the nonlinear anelog of the linear bourcary bus impedance model is 


((PEQI - jQEQ1) + (PB - jQ3)) + (PEQX - jQEQx) 


- {25 ][yeor + Yes + Yeox] ES 


the equations for current injections 


ati Sai) Fai 
TL Ys" YBB| | BB 

YII YB EI 
ox YIB' (YBE+YE9X) | | EB 


(YTL + YEOX)EB - IEQX 


(TEQI + 2) + IEQK = (YEQI + YB3 + YEQK)EB 


| 
| 


Ld 
4H 
td 


< 


bs 


EB 


Le 


a 


(C.33) 


(C.34) 


ue) 


225) 


ooo) 
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EPPENDIX D. MAXIMUM LIFELIHCOD IDENTIFICATION EQUATIOINS 


Given a sequence of input measurerents u(1), u(2),..., u(QX), 
a sequence of corresponding output measurements 2(1), 2(2),..., 2(N), 
a set of confidence coefficients c¢(1), c(2),..., ¢c(N) which indicate 
the relative confidence associated with each input/output set, and 


given the input/output relation 
z2(n) = H u(n) + v(n) (D.1) 


where v(n) is a noise vector having a Gaussian probability distribution 
e e ak e e e e e 
with zero mean and covariance ea) R, it is desired to find the maximum 
likelihood estimates of H and R subject to the following conditions: 
1. The elements of H are independent of one another. 


(e.g. in general H f# H') 
a E $v(n)§ = 0 


3. E $v(n)v' (m)$ -0O fornAm 


ut 


2 Bf u(n)yt(n)f a(n) = 
E fu(n) ut(n)§ 


Wn 
i 
O 


The probability distribution function for the Kx1 vector v(n) is 


Sv! (n)e(n)R7 y(n) 


o[v(n)}] = [(ar)* e7*(n)izl] ~ e (D.2) 


NIH 


NO 


The likelihood function for z(n) and u(n) given H and R where z(n) and 


u(n) are Kxl vectors, and H and R are KxK matrices is 


i + 1 
pfz(n),u(n): BB) = [lar)® o“Rnyint] 2 eFlatr) wR) m3) 


where 





2O3= 


Ifz(n),u(n) 4,2] = 2[2(n) - H u(nlJe(n)R*[2(n) - Eu(n)y 0.4) 


Using trfat to mean the trace of matrix A and using the identity 


a'a = trfa att where a is a vector equation (D.4) can also be written 
- t 
J[z({n),u(n):H,R] = = tr fe(n)a "2(n) - EH u(nJj[2(n) - # a(ny}"} Sy) 


Iet the sets Z, U, and V be such that Z f2(1), (2) 2(r)¢ 
Y= 4u(2), u(2),..., u(t, and Y= {v(2), v(2),..., (ME. 


Because v(n) is assumed to be uncorrelated in time (fv (ny! (m)f =a0 


for n #m) and because v(n) is Gaussian, the joint probability distri- 
bution function for V¥ is the product of the distribution functions for 


MEN (2), 6.04 Ve 


p(v¥) = pfv(1)] pfy(2)}... pfv(y)] 
. X A 4S ve(me(nye v(m) 
fon™ar tetol? 7 
n=] 


The likelihcod function of Z ard U given E and R is then 


p(2,uR) = \Cor)™\at™ Tf ov Xny| * (0.7) 
n= 
where 
ae - : =A \ ' ; 
r(z,u#R) = 4 > teSc(me fem) - £ uln))[z(n) - £ u(n)] $0.8) 
n=] 


It is desired to find the values of H and R which maximize the likelihood 
function p(Z,U:H,R). These values of H and R are the maximum likelihood 


A A 
estimates H and R. If the logarithm of the likelihood function is maxi- 


mized, then so will the function. The log likelihood function is 





OO 


N 
én{p(Z,U:E,R)] — 5 K £n(21) + 5 zs, Keen c (n) 
n=] 
- 5 NAnfRI - 3(Z,U:H,R) (D.9) 


A new function can be cefined 


f(Z,U:H,R) 


N 
-2£n[p(Z,U:H,R)) - NK fn(27) + S K An c(n) 
n=l 


N Al RI + 25(Z,U:H,R) (D.10) 


f£(Z,U:H,R) only has terms which involve H or R, hence, minimizing 

f(Z,U:H,R) is equivalent to maximizing the likelihood function. 
To perform the necessary mathematical manipulations certain 

matrix gradient identities are used from Athans and Schweppe [1 J. 


The matrix gredient for a scalar function of a KxM matrix X is defined 


for use here as 














7a) af (X) 
AX Akay 
i [e(x)] = |----------- (0.11) 
df (x) af (X) 
Ika IX, 


Two properties of trace of a matrix are 
tr fa} = tr far (D712) 


tr fa By = tryB Ay 


The identities uSed are 


J 2 
oe tr fAXt =A (D.14) 
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god, ae 2 ' ryt at 

Sy tr TAX BMG SRK ALEXA (0.15) 
a 3A a BS Se Ep es Gee) 
ox si 5 2 2264 a 
2 eT 

Dx $n{X[ =X (Diyas) 


When £(Z,U,H,R) is at its minimum value, the matrix gradients of 
f(Z,U:H,R) with respect to H and R will be zero. Then the maximum like- 


lihood estimates H end R willimbesthe*solutcdonse fore heendenktor 


a 

OH £(Z, UeheR) = Q (B.18) 

o Y 
and = eR) =O (D.19) 
for 

£(Z,U:H,R) = N 4nfRl + 23(Z,U:H,R) e220), 
and 


Cy 


N t 
(2,U:H,R) = 5 2, tr $¢(n)R7*[ z(n) -H u(n)l a(n) - # u(n)] 01) 
Che 


First, the gradient with respect to H is 


Bae 


a 
<q P(Z,UH WR) = Sy 20(Z,U:H,R) = 0 


z ! 
Fy = te [olde Lal) - Ewe) - # u(n))'t 


N 


~e(n)R7TH uln)zt(n) + ¢(r)R7H w(n)ut (n't (ree) 





=O i= 


Using equations (D.12) and (0.13) 
tr fo(n)R'z(n)ut (ner = tr fe (nz w(n)z" (n)R7g 

= tr $o(n)u(n)z' (na 7H 
tr Jo(mR ME u(n)at(n)$ = tr $o(n)u(n)z! (mR? 


Then equation (D.22) is 


Dy tr } -20(n)u(n)z! (n)R'H + acine a u(n)ut (nz 


and using equetions (D.14) and (D.15) 


N 


x | -2c(n)u(n)z'(n)R7+ + 2c (n)u(n)u! (n)HtB+Y 
n= 


Q 


N N ay 
2\- >, e(n)u(n)z'(n) + >; e(n)u(n)u' (n)E A R 
n=l] 


ni 


-] e e e o a e 
Then unless R~ = QO (which would mean infinite noise covariance) 


N N 
> e(n)u(n)z'(n) = > ec(n)u(n)u'(n)H! 
n=] n=] 


N N 
H 2, e(n)u(n)u'(n) = >) e(n)g(n)u'!(n) 


n= 


and the value of H which maximizes p(Z,U:H,R) is 


P N N a7 
H = |Z o(n)z (mda (n)| | > e(n)u(n)u* (n)| 


N=) Naw 


Taking the matrix gradient of f(Z,U:H,R) with respect to R 


(D 


(D 


(D 


(on 


(D. 


(D. 


23) 


oe) 


25) 


.26) 


27) 


28) 
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3 
oR £(2,U:H,R) = 0 
N 
-N-Se engl + Sa a $c(n)e fz(n) - Hula) [z(n) - # a(n)" 
(Dee0) 
Using equations (D.16) and (D.17) 
N ’ 
Q=ig7 - >. c(n)R ‘f2(n) -Hu(n)}[2(n) - u(n)|' a7 
n=] 
zt yh N : i 
= IR f -% 2 o(n)[z(n) - H u(n)}[z() - 8 atoy't Ro (D. 31) 
ri 


Again unless oa = Q, using equation (2.29) the value of R which mini- 


A N A A 1 
hoe 2 o(n){z(n) - 2 ulnj}[2(n) - £ a(n)] (D.32) 
n= 


A A 
Hence H ard R or eauations (D.29) and (D.32) are the maximum likelihood 


solutions to the problen. 
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APPENDIX E. COMPUTER SIMULATION PROGRAM 
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